Analytical investigation of the exact groundwater divide between rivers beyond the Dupuit–Forchheimer approximation

The groundwater divide is a key feature of river basins and significantly influenced by subsurface hydrological processes. For an unconfined aquifer between two parallel rivers or ditches, it has long been defined as the top of the water table based on the Dupuit–Forchheimer approximation. However, the exact groundwater divide is subject to the interface between two local flow systems transporting groundwater to rivers from the infiltration recharge. This study contributes a new analytical model for two-dimensional groundwater flow between rivers of different water levels. The flownet is delineated in the model to identify groundwater flow systems and the exact groundwater divide. Formulas with two dimensionless parameters are derived to determine the distributed hydraulic head, the top of the water table and the groundwater divide. The locations of the groundwater divide and the top of the water table are not the same. The distance between them in horizontal can reach up to 8.9% of the distance between rivers. Numerical verifications indicate that simplifications in the analytical model do not significantly cause misestimates in the location of the groundwater divide. In contrast, the Dupuit–Forchheimer approximation yields an incorrect water table shape. The new analytical model is applied to investigate groundwater divides in the Loess Plateau, China, with a Monte Carlo simulation process taking into account the uncertainties in the parameters.     

The Bear Brook Watershed in Maine: Multi-decadal whole-watershed experimental acidification

The Bear Brook Watershed in Maine (BBWM) is a long-term research site established to study the response of forest ecosystem function to environmental disturbances of chronic acidic deposition and ecosystem nitrogen enrichment. Starting in 1989, the West Bear (treated) watershed received bimonthly applications of ammonium sulfate [(NH₄)₂SO₄] fertilizer from above the canopy, whereas East Bear (reference) received ambient deposition. The treatments were stopped in 2016, marking the beginning of the recovery phase. Research at the site has focused on soils, streams, and vegetation. Here, we describe data collected over three decades at the BBWM—input and stream output nutrient fluxes, quantitative soil pits and soil chemistry, and soil temperature and moisture.     

2007-2017年中国沙尘气溶胶的三维分布特征及输送过程

利用最新发布的CALIPSO产品，构建了2007-2017年中国沙尘气溶胶的三维分布，并结合HYSPLIT-4模式和再分析数据，探讨了沙尘的三维输送过程。结果表明：中国的沙尘排放源区主要是塔克拉玛干沙漠和巴丹吉林沙漠，沙尘气溶胶出现频率分别为60%和35%。塔克拉玛干沙漠排放的沙尘主要（50%~70%）停留在源地0~6 000 m高度，少部分向东输送至甘肃和内蒙古；巴丹吉林沙漠排放的沙尘则主要向东输送。中国沙尘排放量在春季最大，向东输送最强；夏季，东亚夏季风限制了沙尘向东输送；秋季，沙尘排放减弱，输送强度和夏季相当；沙尘排放量在冬季最小，输送最弱。夏季，沙尘在输送过程中可被抬升至高度5 000 m以上，春季次之，秋、冬季的沙尘主要在低层大气输送。沙尘在向东输送的过程中被抬升并和当地人为污染物混合变为污染性沙尘，华北地区污染性沙尘出现频率高达30%；输送到海洋的沙尘也会与洋面上（0~3 000 m高度）的海盐气溶胶混合，出现频率约为10%。     

Seasonal connections between meteoric water and streamflow generation along a mountain headwater stream

In snowmelt-driven mountain watersheds, the hydrologic connectivity between meteoric waters and stream flow generation varies strongly with the season, reflecting variable connection to soil and groundwater storage within the watershed. This variable connectivity regulates how streamflow generation mechanisms transform the seasonal and elevational variation in oxygen and hydrogen isotopic composition (δ¹⁸O and δD) of meteoric precipitation. Thus, water isotopes in stream flow can signal immediate connectivity or more prolonged mixing, especially in high-relief mountainous catchments. We characterized δ¹⁸O and δD values in stream water along an elevational gradient in a mountain headwater catchment in southwestern Montana. Stream water isotopic compositions related most strongly to elevation between February and March, exhibiting higher δ¹⁸O and δD values with decreasing elevation. These elevational isotopic lapse rates likely reflect increased connection between stream flow and proximal snow-derived water sources heavily subject to elevational isotopic effects. These patterns disappeared during summer sampling, when consistently lower δ¹⁸O and δD values of stream water reflected contributions from snowmelt or colder rainfall, despite much higher δ¹⁸O and δD values expected in warmer seasonal rainfall. The consistently low isotopic values and absence of a trend with elevation during summer suggest lower connectivity between summer precipitation and stream flow generation as a consequence of drier soils and greater transpiration. As further evidence of intermittent seasonal connectivity between the stream and adjacent groundwaters, we observed a late-winter flush of nitrate into the stream at higher elevations, consistent with increased connection to accumulating mineralized nitrogen in riparian wetlands. This pattern was distinct from mid-summer patterns of nitrate loading at lower elevations that suggested heightened human recreational activity along the stream corridor. These observations provide insights linking stream flow generation and seasonal water storage in high elevation mountainous watersheds. Greater understanding of the connections between surface water, soil water and groundwater in these environments will help predict how the quality and quantity of mountain runoff will respond to changing climate and allow better informed water management decisions.     

Variations in glacier volume and snow cover and their impact on lake storage in the Paiku Co Basin,in the Central Himalayas

Glaciers and snow cover are important constituents of the surface of the Tibetan Plateau. The responses of these phenomena to global environmental changes are sensitive, rapid and intensive due to the high altitudes and arid cold climate of the Tibetan Plateau. Based on multisource remote sensing data, including Landsat images, MOD10A2 snow product, ICESat, Cryosat-2 altimetry data and long-term ground climate observations, we analysed the dynamic changes of glaciers, snow melting and lake in the Paiku Co basin using extraction methods for glaciers and lake, the degree-day model and the ice and lake volume method. The interaction among the climate, ice-snow and the hydrological elements in Paiku Co is revealed. From 2000 to 2018, the basin tended to be drier, and rainfall decreased at a rate of −3.07 mm/a. The seasonal temperature difference in the basin increased, the maximum temperature increased at a rate of 0.02°C/a and the minimum temperature decreased at a rate of −0.06°C/a, which accelerated the melting from glaciers and snow at rates of 0.55 × 10⁷ m³/a and 0.29 × 10⁷ m³/a, respectively. The rate of contribution to the lake from rainfall, snow and glacier melted water was 55.6, 27.7 and 16.7%, respectively. In the past 18 years, the warmer and drier climate has caused the lake to shrink. The water level of the lake continued to decline at a rate of −0.02 m/a, and the lake water volume decreased by 4.85 × 10⁸ m³ at a rate of −0.27 × 10⁸ m³/a from 2000 to 2018. This evaluation is important for understanding how the snow and ice melting in the central Himalayas affect the regional water cycle.     

Benthic foraminiferal fauna and sediment provenance of Eocene syn-rift sequences in Taiwan: implication for onset of Asian epi-continental marginal seas off China coast

Marine Geophysical Research - Early Cenozoic rift basins developed commonly on the Mesozoic basement along the SE Asia Continent. However, Eocene–Oligocene sequences were only exposed widely...     

Creating an isotopically similar Earth–Moon system with correct angular momentum from a giant impact

The giant impact hypothesis is the dominant theory explaining the formation of our Moon. However, the inability to produce an isotopically similar Earth–Moon system with correct angular momentum has cast a shadow on its validity. Computer-generated impacts have been successful in producing virtual systems that possess many of the observed physical properties. However, addressing the isotopic similarities between the Earth and Moon coupled with correct angular momentum has proven to be challenging. Equilibration and evection resonance have been proposed as means of reconciling the models. In the summer of 2013, the Royal Society called a meeting solely to discuss the formation of the Moon. In this meeting, evection resonance and equilibration were both questioned as viable means of removing the deficiencies from giant impact models. The main concerns were that models were multi-staged and too complex. We present here initial impact conditions that produce an isotopically similar Earth–Moon system with correct angular momentum. This is done in a single-staged simulation. The initial parameters are straightforward and the results evolve solely from the impact. This was accomplished by colliding two roughly half-Earth-sized impactors, rotating in approximately the same plane in a high-energy, off-centered impact, where both impactors spin into the collision.     

短脉冲激光在实验室天体物理方面的研究进展

随着啁啾脉冲放大技术(Chirped Pulse Amplification,CPA)的飞速发展,激光功率密度实现了飞跃式的提升,利用短脉冲激光开展实验室天体物理研究的条件日趋成熟。短脉冲激光与靶相互作用可以产生相对论粒子(正负电子、质子、中子等)和高能电磁辐射(X射线、γ射线),这些粒子和辐射的产生过程与天体中的某些物理现象相对应。介绍了利用短脉冲激光开展实验室天体物理的研究进展,重点介绍了核天体物理、空心离子光谱学、相对论无碰撞冲击波以及等离子体状态诊断等。     

Three-dimensional finite element simulation of fracture propagation in rock specimens with pre-existing fissure(s) under compression and their strength analysis

A FORTRAN program, consistent with the commercially available finite element (FE) code ABAQUS, is developed based on a three-dimensional (3D) linear elastic brittle damage constitutive model with two damage criteria. To consider the heterogeneity of rock, the developed FORTRAN program is used to set the stiffness and strength properties of each element of the FE model following a Weibull distribution function. The reliability of the program is assessed against available experimental results for granite cylindrical specimens with a throughgoing, flat and inclined fissure. The calibration procedure of the material parameters is explained in detail, and it is shown that the compressive to tensile strength ratio can have a substantial influence on the failure response of the specimens. Numerical simulations are conducted for models with different levels of heterogeneity. The results show a smaller load bearing capacity for models with less homogeneity, representing gradual coalescence of fully damaged elements forming throughout the models during loading. The maximum load bearing capacity is studied for various combinations of inclination angles of two centrally aligned, throughgoing and flat fissures of equal length embedded in cylindrical models under uniaxial and multiaxial loading conditions. The key role of the compressive to tensile strength ratio is highlighted by repeating certain simulations with a lower compressive to tensile strength ratio. It is proven that the peak loads of the rock models with sufficiently small compressive to tensile strength ratios containing two throughgoing fissures of equal length are similar, provided that the minimum inclination angles of the models are the same. The results are presented and discussed with respect to the existing experimental findings in the literature, suggesting that the numerical model applied in this study can provide useful insight into the failure behaviour of rock-like materials.